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We formulate a momentum-shell renormalization group (RG) procedure that can be used in theo- 
ries containing both bosons and fermions with a Fermi surface. We focus on boson-fermion couplings 
that are nearly forward-scattering, i.e. involving small momentum transfer (g ~ 0) for the fermions. 
Special consideration is given to phase space constraints that result from the conservation of mo- 
mentum and the imposition of ultraviolet cutoffs. For problems where the energy and momentum 
scale similarly (dynamic exponent z = 1), we show that more than one formalism can be used and 
they give equivalent results. When the energy and momentum must scale differently (z 7^ 1), the 
procedures available are more limited but a consistent RG scheme can still be formulated. The 
approach is applicable to a variety of problems, such as itinerant-electron magnets and gauge fields 
interacting with fermions. 
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I. INTRODUCTION 

Although the theory of seahng and renormalization 
has profoundly affected our conceptual understanding of 
many-body systems, its calculational framework is im- 
perfect and continually evolving. In the end, we are in- 
terested in how couplings flow under changes of scale, 
but a variety of distinct procedures exist, each with its 
own advantages and drawbacks. An incomplete list of 
the assortment of programs includes the multiplicative 
RG, real space decimation, functional RG, exact RG, 
flow equations, and various flavors of e-expansion, such 
as the classic minimal subtraction which expands around 
d = 4, or expansions around some other parameter, such 
as the deviation of the range of the interaction from a 
suitable reference value. Each method has its own lim- 
its of practicality, ease of use, and range of problems to 
which it may be usefully employed. Wilson's momentum- 
shell approachii^ is an especially popular method in the 
context of condensed matter problems. However, in the 
early 1990's a few people recognizedM;^ that the stan- 
dard momentum-shell procedure must be modified for 
problems involving a Fermi surface. A campaign soon 
followed attempting to understand Fermi liquid theory 
from an RG perspective. An excellent and influential 
summary of the pure fermion RG can be found in-. 

Another indication that the RG for fermions required 
more scrutiny came from the study of quantum criticality 
in itinerant electron magnets. The usual Hertzian ap- 
proach^ uses an auxiliary (Hubbard-Stratonovich) field 
to decouple the fermion-fermion interaction, thus allow- 
ing fermions to be completely integrated out. The re- 
sulting effective theory is then expressed in terms of 
the remaining bosonic auxiliary field, to which standard 
bosonic RG techniques can be employed. However, be- 
cause the fermions are gapless, the process of integrat- 
ing them out may introduce non-analyticities in the cou- 
plings among the remaining bosonic modesii^. It would 
therefore be important to devise an RG scheme capa- 
ble of simultaneously handling both bosons and fermions 
with a Fermi surface. 



Besides the critical itinerant magnets, a mixed 
fermionic-bosonic RG formalism would be quite useful 
for an assortment of problems. For example, in the con- 
text of a gauge field coupled to fermions, several au- 
thors^iifliiiiiH have developed their own RG schemes for 
counting dimensions in these mixed theories. All have in 
common the subdivision of the Fermi surface into a large 
number of patches, but results vary and despite the inter- 
vening 15 years since this pioneering work, little progress 
has been made. The importance of the gauge-fermion 
problem is historically linked to an interesting path to 
non-Fermi liquid behavior-i^^i^. More recently, effective 
gauge theories have appeared in a number of additional 
contexts in condensed matter physics^^. 

We should mention in passing a growing body of work 
on the functional RG which may be adapted for mixed 
boson-fermion theorie a^^i^^'^^i^^'^°i^^ . This typically re- 
quires blending with computational methods and may 
prove to be a useful framework for understanding realis- 
tic material bandstructures. Our aim here is rather more 
modest, which is to develop an RG scheme for mixed the- 
ories with a high score in the "ease of use" category. This 
was the chief virtue of the original Wilsonian RG which 
could quickly identify the relevant and irrelevant opera- 
tors with minimal effort. One emphasis of this paper will 
be to carefully consider how to extend Shankar's scheme^ 
to include bosons while maintaining the easy-to-use spirit 
of the Wilsonian approach. 

Our primary motivation to consider these issues came 
from the context of magnetically ordered phases of some 
itinerant systems. In the case of an antiferromagnetic 
state of a Kondo lattice, the bosonic magnons, described 
by a quantum non-linear sigma model, are coupled to 
the fermionic quasiparticles near a Fermi surface. In this 
problem, energy and momentum scale the same way; the 
dynamic exponent z = 1. The ferromagnetic counter- 
part features z = 3. The RG analysis plays an essential 
role in understanding the Fermi surfaces in these systems 
and has been briefly described in our earlier works22i^i2i. 
The purpose of the present work is to explain the details 
of the method in considerable detail with the hope that 
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the method will be adapted to problems in new physical 
contexts. 

The remainder of the paper is organized as follows. In 
section we remind the reader of the essential points 
of the bosonic Wilson-Hertz scaling. Section [Till quickly 
moves on to discuss scaling in fermionic systems, largely 
paraphrasing what has already been done but emphasiz- 
ing a slightly different perspective. The next section. ITVl 
describes a method to properly scale in mixed theories 
when energy and momentum can be given the same scal- 
ing dimension, i.e. when z = 1. This is closest in spirit 
to the Shankar approach, but cannot be generalized to 
z ^ 1. Some of these problems are discussed in sections IVl 
and IVII In section IVTl we present an alternative method 
for arbitrary z that is perhaps less intuitive than that of 
section IIVI but has the advantage of being generalizable 
to z 7^ 1 while at the same time yielding identical results 
when z = 1. 



II. BOSON SCALING 



review those elements crucial to the comparisons we wish 
to make later with the fermionic RG. 

Consider a d-dimensional integral in momentum space 
with a cutoff to high-energy and therefore large-g modes; 
this is denoted by A. Let us separate out a thin shell 
of high energy modes in the range A/s < g < A, where 
s > 1. 







A/s pA 

q'^-^dq+ / q-^-^dq 
J A/s 



Here, g = is the radial coordinate in (hyper)spherical 
coordinates, and d'^^^il^ represents the measure for inte- 
gration over all angular variables in g-space. We have ig- 
nored factors of 27r. Mode elimination amounts to simply 
throwing away the shell integral. To regain the original 
form of the action only a trivial rescaling of the radial 
coordinate is needed: 



The problem we are concerned with can be decomposed 
into bosonic, fermionic, and interaction terms: 



+ 5^ 



(1) 



The bosonic and fermionic pieces can be further divided 
into quadratic and quartic pieces. 



5" 



(2) 
(3) 



Theories based upon or alone have already been 
subjected to momentum-shell RG analyses; see, for 
example^i^ and^. In this section, we review the Wilson- 
Hertz scaling procedure for bosons, so we are only con- 
cerned with S^. 

In the most general case, the quadratic part of the 
action can take several different forms depending on the 
value of z. For example. 
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The bosons might represent acoustic phonons, magnons, 
photons, or some collective mode of an underlying 
fermionic theory that results after "integrating out" the 
fermions with an auxiliary field. At this point, we need 
not be specific. All that matters is that we must design 
the RG scheme in such a way that ^2 remains invariant. 
The Wilsonian RG for bosons is well-knowiii so we only 



q' = sq 



(5) 



This defines the scaling dimension of momentum. In the 
customary notation, we use square brackets to denote the 
scaling dimension of any quantity according to 



A' 



(6) 



where A' is measured in units s''^' times smaller than 
the units of A. We call [A] the scaling dimension of A. 
This notation differs from another frequent convention 
which may claim, for instance, [volume] = where L is 
some length scale. We prefer our notation since it means 
a coupling g is relevant when [g] > 0, irrelevant when 
[g] < 0, and marginal when [g] =0. 

In this notation, an equivalent statement to equation 
([5]) is simply 



(7) 



Using this form of momentum scaling in the integral leads 
to 

J d'-'n,^j\-<-^-'^q'^-h-'dq' = s-" {d^q)' 

We conclude that the scaling dimension of the measure 
is given by 

[dS] = d[q] = d (8) 
Note that rescaling the radial variable, q, is the same 



as rescaling all the components of q since q ^ y la ■ 
For this to be consistent with q' = sq, we must have 
q'^ — sqa for all components a G {x, ?/, z, ...}. This is 
an important difference from the fermionic case to be 
discussed later and results from the simple fact that the 
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coordinate origin here is a single point rather than an 
extended surface. 

Let us apply this mode elimination and rcscaling to 
the quadratic part of the boson action, taking the case 
z = 3 as an example. 



'd[qh 



d'^q'duj'\(f>{s~^l^q',s-^''kLo')\'^ 



q' 



(9) 



where we discarded high energy modes and used equa- 
tion ^ and Lj' = st'^lw. We wish to scale the terms 
in the parentheses identically, so we must have ~'2[q\ — 
— [lu] + [q] thus fixing the relationship between the scaling 
dimensions of energy and momentum: 



[q] = H/z 



(10) 



where our example considered z = 3 explicitly. 

The final step is wavefunction renormalization which 
can be implemented by defining a new field (p' according 
to: 

iLu') = .s-{(''+2)M + M}/2^(,s-Mq', s-Mjc^') (H) 

Equivalently, we might say that the boson field has a 
scaling dimension given by 



[0] = -M(d + z + 2) 



(12) 



where we used [uj] ~ z[q]. Although customary, equa- 
tion p2|) is slightly misleading. The replacement of 
by s~['^l0' in analyzing interaction terms should only be 
done when the arguments of the field transform accord- 
ing to equation (jlip . Equation (|6ip provides an example 
where the arguments of the field are transformed in a 
very different fashion. From equation (fTTj) we see that 
the boson field appears to take the form of a generalized 
homogeneous function. We do not delve into this issue 
further, but merely note that equation (|lip is a very spe- 
cific type of substitution that needs to be implemented 
in this strict form. Scale invariance of S2 has imposed 
a transformation property on the field, specified in (|lip . 
under the particular coordinate transformation q' ~ s^'^q 
and w' = fil'^lw with [g] = [a;]/z. 

Now that we know how to scale momentum from equa- 
tion ([7]), energy from equation (fTO|) . and the field from 
equation ([TT|) . we are ready to analyze the four-boson 
interaction term: 
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d"^ q4^dLJ4^d'^ q^duj^d'^ q2diL!2d'^ qiduJi 



X(t){q4,iuJ4)(t){ci3, iu)3)(l){q2 , iuj2)4>iqi , iuJi) 

xe(A - |q4|)e(A - |<f3l)e(A - |<f2|)e(A - m 

x5'-'^\qi + (fa - 92 - qi)S{uj4: + UJ3 - UJ2 - i^i) 

(13) 



The ^-functions enforce the conservation of energy and 
momentum while the 0-functions define the cutoffs for 
the effective field theory (in principle, energy cutoffs 
should also be written, but this is understood). 

To determine the scaling dimension of ut at the tree 
level, we first separate the integrations into low and 
high energy modes [i.e. 9(A — \qi\) = 6(A/s — \qi\) + 
Q{\qi\ — A/s)8(A — \qi\)], then discard the high energy 
shell. There is some freedom in choosing the shape of the 
shell which can take some curious forms for the purpose 
of simplifying calculations. See the discussion by Herta^. 

After rescaling according to ([7]), pUj) . and ([H]), we find 



A~d-z 



Ub 



d'^q'^duj'sd'^q'^duj'^d'^q'^duj'i 



X(l){q'4^,iu'^)(p{q'^,iuj!^)(j){cf2,ii^'2)<P{q'i,'i(^'i) 

xe(A- |9i|)e(A- |g'3i)e(A- |g-|)e(A- kii) 

XS^"^ (ql + (fg - (f^ H +^3-^2- ^'1) 

(14) 



which tells us that = ^Ub, or equivalently 

[ub] = A-{d + z) 



(15) 



This yields a quick way to determine when the four- 
boson interaction term Ub J 0^ is relevant or irrelevant 
based on the dimensionality of the problem and the value 
of z. Historically, this result provided some early intu- 
ition about quantum phase transitions which can behave 
like classical phase transitions but in a different num- 
ber of effective dimensions: dcs ~ d + z. Although the 
theory was originally devised to address questions about 
itinerant quantum critical magnet a^'^^i^^ , some problems 
have been encountered with this approach^i^. Part of the 
problem could be that the theory is completely bosonic, 
despite the underlying fermionic nature of the system. It 
is therefore desirable to develop an RG formalism that 
includes fermions with a Fermi surface. 



III. FERMION SCALING: SHANKAR'S RG 



by 



For fermions, the quadratic part of the action is given 



si = I d^Kde ip {ie - iP 



(16) 



To define a scaling scheme that leaves scale invari- 
ant, we now review the formulation of the fermionic RG^. 
We shall use Shankar's notation and label momenta mea- 
sured with respect to the Brillouin zone center with a cap- 
ital letter K — {Kx,Ky, ...). In contrast to the bosonic 
case, low energy modes live near an extended surface (the 
Fermi Surface) rather than a single point (the Brillouin 
Zone center). For a spherical Fermi surface, a high energy 
cutoff can be implemented on i?-integrals as follows: 



d^K 



Kp+A 



K'^-^dK 



Kf-A 



4 



where A is an ultraviolet cutoff, but we still insist A <^ 
Kp. Here, represents the measure for inte- 

gration over all angular coordinates in J?-space, while 
K = \K\ is the radial coordinate. Usually, we work at 
fixed fermion density which, by Luttinger's theorem, dic- 
tates that we design our scaling scheme in such a way 
that the Fermi volume remains invariant. To preserve 
the Fermi surface under rescaling we cannot simply scale 
the radial coordinate as we did in the bosonic case. To see 
this, observe that after mode elimination the expression 
we wish to rescalc is given by 



A/s 



d'^-^ng. / K'^-^dK (17) 



Clearly, no simple rescaling of K will return the integral 
to its original form. This is the principle disparity be- 
tween the fermionic and bosonic RG. To make progress 
we define the lower case letter k = \K\ — Kp. Note 
that fc = corresponds to = — /i = since 

~ ^ 2m'' ^ "-f(^ ~ Kp) = vpk. Small k corre- 
sponds to low energy whereas small K does not. Such a 
change of variables greatly facilitates rescaling. 



d'^-^ng 



A/s 



{Kp + kf-^dk 



-A/s 



Kp^ I d<'-^^g 



Kp' I d'^-'ng 
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-A/s 
A/s 

-A/s 
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dk 



dk 



(18) 



We have neglected certain terms above for two reasons: 
they are of order A/Kp relative to what has been kept, 
and they are less relevant in the RG sense. To see the lat- 
ter, note that the integral can be restored to its original 
form with the simple rescaling k' — sk. This determines 
the scaling dimension 



[k] = 1 



(19) 



Note that the variable k is not a vector, nor is it a radial 
coordinate since it can take negative values. Later, we 
will discuss another scheme, which we call "patching," 
that decomposes the momenta into components parallel 
(fc||) and perpendicular {k±) to the Fermi surface nor- 
mal. To make later contrast with the patching scheme of 
section IVTl which uses local coordinates for each patch, 
we will call the present approach the "global coordinate" 
scheme. 

To further emphasize the dissimilarity between the 
fermionic and bosonic cases, observe that after the rescal- 
ing of equation 



A/.S 



d'^K 



A 

Kp' I d^-'ilg I 



'dk' 



A 



{d'^K)' 



(20) 



which implies that, effectively, 

[d'^K] = 1 



(21) 



This stands in sharp contrast to the bosonic case in equa- 
tion ([8|). Here, the angular variables arc truly untouched 
after rescaling which is necessary to maintain the Fermi 
surface. Unfortunately, the straightforward transforma- 
tion k' = sk does not translate into a simple transforma- 
tion on the components of K. Care must therefore be 
exercised to write all expressions in terms of k before the 
scaling procedure can begin. For example, after mode 
elimination and rescaling of energy and momentum, the 
quadratic part of the fermionic action is given by: 

Sl (X s-^ j dk'd€'4>{Kp + s^'k' ,s~'ie') 

X [s-^ie' - WFS"[''lfc'j iP{Kp + s-'k', s'^ie') 

If we wish to scale both of the terms inside the square 
brackets identically, we must choose 



[k] - H 



(22) 



thus fixing the relationship between the scaling dimen- 
sions of energy and momentum. For convenience we can 
set this value equal to 1, as in equation (|19[) . Compare 
this to equation (fTO|l . 

In order to make invariant to the RG transforma- 
tion we must demand that the fermion field obeys: 

s-'^^^i^{Kp + s-'k',s-'ie') = i''{Kp + k',ie'){23) 

where we have not explicitly written the dependence of 
■0 on angular variables since these do not scale. Equation 
(|23p tells us two important things. First, the dimension 
of the fermion field is simply: 



-3/2 



(24) 



Second, the RG transformation of the fermion field does 
not take the form of a generalized homogeneous function 
as was the case for the bosonic field; see equation (|lip . 
The momentum argument of the fermion field K has a 
magnitude equal to the Fermi wavevector plus a small de- 
viation: K = Kp + k. Only the deviation k scales, while 
Kp remains constant. This important difference from 
the bosonic case will be discussed further in section [V] 

The story so far seems relatively elementary, but the 
true subtleties materialize when we try to determine the 
dimension of the -0^ coupling function Uf based on the 
dimension assignments required to make scale invari- 
ant. The quartic part of the action can be written^ 

5/ = n y '^''^'^ J de~5^''\Ki +K2-Ks~ K^) 



X5(ei -I- £2 - £3 - £4) 

x0(4)V;(3)0(2)V(l)^^/(4,3,2,l) 



(25) 
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The 5-functions explicitly enforce the conservation of en- 
ergy and momentum (up to a reciprocal lattice vector). 
We might integrate one of the energies and momenta, 
say (1^4,64), against the delta function to yield an inte- 
gral over three independent sets {Ki,ei), {K2,^2), and 
(i?3,e3). 



Si = 



cfiKi / (ie,V(l + 2 - 3)^/'(3)V'(2)V'(l) 



n 

i=l 

xu/(l + 2 - 3,3,2, 1) (wrong) 



(26) 



But this expression is not quite right. The problem is 
that not all momentum-conserving processes should be 
included in the low-energy effective field theory. We must 
respect the cutoff imposed on the quadratic part of the 
action, which only allows excursion into states within 
a distance ±A of the Fermi surface. Imposing a cutoff 
amounts to constraining the momentum integrals. Un- 
til now, we have implemented the cutoff constraints by 
writing them explicitly in the limits of integration, but 
let us re-express them as 



(27) 



where, as usual, ki = \Ki\ — Kp- With all momentum 
integrals written in this way, we can safely use the 5- 
functions to eliminate one variable, say and £4. 

5/ = n /" d'-K^ I de.,i,{l + 2 - 3)Vi(3)V'(2)V'(l) 

i=i 

xw/(l + 2-3,3,2,l) 

xe(A- |fci|)e(A- |fc2|)e(A- |fc3l)e(A- 1/C4I) 

= n / '^''^^ / '^^^^^(l + 2 - 3)^(3)V'(2)V(l) 



(28) 



xw/(l + 2-3,3,2,l)e(A- I/C4I) 



The constraints on Ki, K2, and K3 have been put back 
in the limits of integration, but we have the additional 
constraint I/C4I < A where 



IC4 



(29) 



Once we have conserved momentum, is no longer an 
independent variable, so we use the notation /C4 to rep- 
resent the combination of variables specified in equation 

We can implement the constraint embodied in 0(A — 
I/C4I) in a number of ways. One way is to allow A'l 
and K2 to range anywhere inside the annuli defined by 
—A < ki,k2 < A, but restrict as appropriate to sat- 
isfy I/C4I < A. The outcome of a proper phase space 
analysis shows that once Ki and K2 have been chosen, 
the angle for K3 is highly constrained^. 



To see this in more detail, observe that to leading order 
in A/A>, 



^f(|A|-1) 



(30) 



where A = A'l + K2 — K3, and where the Ki are unit 
vectors, each pointing in the direction of Ki. Note that 
A is not itself a unit vector since 



lAI = V2 



- + Ki- K2 - Ki - K3 - K2 ■ K3 



1/2 



(31) 



a result we will use in section [V] After mode elimination, 
the momentum integrals become: 



3 » 

Y[ / d'^KS (a/.5 - Kf 

i—l 

3 » 



lAI -1 



lAI-l 



(32) 



Simply rescaling k[ = ski is not sufficient to regain the 
original form of the action for generic values of ki. The 
obvious snag is the annoying way the 0-function trans- 
forms. For general values of the momenta ki, the O- 
function is clearly not invariant to the renormalization 
group transformation. Consequently, we are not techni- 
cally entitled to compare the coupling before and after, 
so we do not know the RG flow. The way out of this 
dilemma is first to understand the circumstances under 
which the 0-function is invariant, and then to see what 
might be happening for more generic cases by considering 
a soft cutoff. 

First, note that when |A| = 1 (i.e. /C4 = 0) the O- 
fimction is always form- invariant since 0(A) = 0(A/s). 
The condition |A| = 1 can be fulfilled in three different 
ways: 

(i) Ki = K3 and K2 = i?4 (33) 

(ii) K2 ^ K3 and Ki = Ki (34) 
(ui) A'l = -K2 and A'a = -A^4 (35) 

For these values of the momenta, the rescaling fc^ = 
ski works flawlessly because the 0-function is form- 
invariant under these restrictions. We are now al- 
lowed to compare the coupling before and after. Since 
[dkidk2dk3deide2de3] = 6 and [■0^] = —6, we conclude 
that, at the tree level, the most relevant pieces of m/ are 
marginal. This important result is at the heart of Fermi 
Liquid Theory, but is expressed by the simple equation: 







(36) 



In Shankar's notation, cases (i) and (ii) correspond to 
Uf — zLF and case (iii) Uf = V. It has also been shown^ 
that case (i) remains marginal beyond the tree level, 
while loop corrections in case (iii) lead to a marginally 
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relevant coupling for certain angular momentum chan- 
nels, indicative of the BCS instability. 

Let us understand in more detail the circumstances 
under which the 0-function is always form-invariant. In 
particular, we want to stress that the the condition IC4 = 
is conceptually different from the limit A/Kp ^ 0. To 
see this, let us rewrite the equation /C4 = as follows: 



device, the rescaled cutoff for arbitrary ki values becomes 



(37) 



Note that equation ((29|) is slightly more accurate than 
([50)1 . Next, define P = Ki + K2 which obviously gives 



\K3-P\ = Kf 



(38) 



This says that the vector joining the tip of to the 
tip of P must have magnitude precisely equal to Kp- 
Figure [T] depicts the situation. Geometrically, the choices 



Ki must stay within the annulus 




- P 



— P must be a point on this circle 

FIG. 1: Once Ki and K2 have been chosen the conservation 
of momentum and the requirement that all Ki respect the 
cutoff of the field theory strongly constrains the phase space 
available to K3. Shown here is the limit /C4 = 0. Even in this 
limit, there is some flexibility in the choices available to K^, 
as depicted by the thick gray lines. While the magnitude of fcs 
can still fall anywhere in the range —A < fca < A, the angles 
available to K^ are highly limited. If we were to further take 
the limit KjKp —> 0, the gray lines would shrink to points. 
Note that P is defined as the sum of Ki and K2, but the 
latter are not drawn to avoid clutter. 

available to A'3 once A'2 and A'l have been selected are 
given by the thick gray lines in the figure. Notice that 
while can still take any values —A < k^ < A within the 
annulus, the angle of K3 has become highly constrained. 
However, it is clear that even when /C4 = the value of 
A/Kp can still be nonzero. 

We now know the dimension [uf] when we restrict to 
IC4 = 0. However, the three cases corresponding to IC4 = 
constitute only a small portion of {Ki, K2, K3)-space. 
To see what happens to the coupling function u(3, 2, 1) 
for other values of momenta. Shankar had the insight to 



e A - sKf 



,-s7Va||A|-1| 

,-Wa||A|-1| -{.,-1)Wa||A|-1| 



where we have defined the large parameter A^a = Kf/ A 
(generally, we have the hierarchy k < A Kp, which 
means A^a ^ 1). We choose to write the cutoff in this 
way because then clearly when |A| = 1, corresponding 
to the three cases listed above, the cutofi^ becomes a 
simple factor of unity. For any |A| 7^ 1, which means 
all other values of the ki, the cutoff ^ in the limit 
A'a — > 00 provided s > 1. While we do not know how 
the coupling function it(3, 2, 1) scales for values of the mo- 
menta where |A| ^ 1. it does not matter because these 
couplings will be exponentially suppressed in the limit 
l/AT^ = A/Kp 0. 

Note that the condition |A| = 1 is simply the state- 
ment that /C4 should not scale. Indeed, it means /C4 = 0. 
Only when |A| = 1 is the relation /C4 = S/C4 satisfied, 
albeit trivially. In fact, that is how we identified the con- 
dition |A| = 1, being the only combination of K^, K2, 
and Ki where 0(A — I/C4I) can be rescaled to take its 
original form after mode elimination. This useful inter- 
pretation will be used again later when we extend the 
formalism to include bosons. 

Before moving on, we need to make another observa- 
tion about the pure fermion RG that will be important 
to later generalizations. We have shown how to find the 
dimension of the coupling function uj(3,2,l) for those 

values of momentum that satisfy IC4 = {i.e. |A| = 1) 
corresponding to forward, exchange, and Cooper scatter- 
ing. To be pedantic, this phase space restriction should 
be incorporated into the form of the coupling: 



f _ 



n 



d'^K, / de^^il + 2 - 3)i^(3)V'(2)i^(l) 



xufil + 2 - 3, 3, 2, l)e(A - \IC4\)S{IC4) (39) 

Note that 0(A) = 1 always since A > 0. As seen in 
figure[l] the insertion of <5(/C4) does not affect the freedom 
of Ki or K2 at all, nor does it affect the magnitude of 
K3 so long as — A < |/?3| — Kp < A. However, the angle 
of K3 is highly restricted to the two gray regions of the 
figure as a resulting of inserting (5(/C4). We may therefore 
implement the constraint (in d ~ 2) by: 



J(/C4) ^ 6i\A\ - 1)/Kp 

= [s{e3-9i) + s{e3-92)]/Kp 



(40) 



employ a soft cutoff: 0(A — \ki\) w e 



-Ife.l/A 



Using this 



A similar expression can be written in d = 3. Since angles 
do not scale in this scheme, whether or not we insert this 
factor into S( will have no effect on the value of the 
dimension of Uf. Shankar's result of marginality, [uf] = 
0, still holds. We mention this issue because generalizing 
the method to include bosons will not result in so happy 
a circumstance. We turn to this case next. 
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IV. BOSON+FERMION SCALING 

We are finally ready to incorporate bosons. Consider 
the following interaction term involving two fermions and 
one boson: 

S'/ = J d''Kid''K2d''qg{Ki,K2,q)^K,^il^^g 

xS^'^^K^ - Ki - '?)e(A - |fci|)e(A - ihi) 
xe(A-lgl) (41) 

g is the coupling function which plays the same role as 
in the 4-fermion problem. For simplicity we have sup- 
pressed frequency integrals and assumed Ah ~ A/ ~ A. 
To conserve momentum wc have two choices: use the d- 
function to eliminate a fermionic momentum Ki, or the 
bosonic momentum q. This gives either 



d'Kd'^q V5^+,^Vj?09 5(^,'?)e(A-|/C2|) (42) 



d''Kid''K2 ['^g^tPji^<l>K,-K, g{K2,K\)e{A- \Q\) 

(43) 

where some of the cutoff constraints have been put back 
in the limits of integrations, and where we have defined 



IC2 = \K + q]-KF 

Q = K1-K2 



(44) 
(45) 



This is analogous to what we did for the pure fermion 
problem; see equations ([28]) and (j29|) . Note that because 
we integrated against the delta functions, momentum and 
energy are already explicitly conserved. In equation 
K2 is no longer an independent variable, so we use the 
symbol IC2 to represent the combination of variables spec- 
ified in equation (|44p . Likewise, qis not an independent 
variable in equation ([43]), so we use Q as shorthand for 
the momentum transfer, as specified in equation (j45|) . 
This mirrors the development of the pure fermion case. 

Unlike the pure fermion problem, we now appear 
to have two different choices for expressing the boson- 
fermion coupling. Equation ()42p involves the boson- 
fermion coupling function g{K,q), while equation (j43p 
contains g{K2,Ki). We defer a discussion of the resolu- 
tion of this choice to section |Vl Here, we simply point 
out that a consistent scheme can only be found for equa- 
tion (im, and we adopt this choice for the remainder of 
this section. 

Although momentum is conserved, just like the pure 
fermion case, not all momentum conserving processes are 
allowed because some might fall outside the high-energy 
cutoffs. We must further restrict the coupling function 
g with the constraint 0(A — |/C2|)- Unfortunately, this 
quantity only scales in a simple way when z ~ I. Let us 
briefly explain the problem. 



Recall from the form of that we have the relation 
[k] = [e] , while S2 demands [q] = [uj]/z for general values 
of z [see equations ([22]) and (fTO|) ]. In addition, since 
we want to scale fermions and bosons at the same time, 
we choose to scale the energies the same way, that is: 
[uj] = [e]. For convenience, we set the scaling dimension 
of energy to unity: [uj] = [e] = 1. Any other value would 
change all scaling dimensions by the same multiplicative 
factor, but their relative dimensions would be unaffected. 
Using this prescription we find 

[e] = [k] = [u;] = l 

[q] = HA-i 
z 



m = - 

[0] - - 



2z 



(46) 



Mode elimination and rescaling according to this scheme 
leads to the following interaction term (we reinstate the 
energy integrals): 



bf _ 



d'^q'dk'd'^-^flj^de'duj'^P'^P'tj)' 
xe{A/s-\IC2\) (47) 



s: 



The reason why we have A/s rather than A/s^/^ is be- 
cause this constraint comes from the restriction on the 
momentum integration of /s2 = |/^2| — Kp in equation 
(|¥T|) . which scales like a fermion. 

Let us rewrite the expression involved in the O- 
function: 



I/C2I = \K + q\-KF 

- [{Kp + kf + q^ 



2{Kp + k)q cos O^q 
-,1/2 



1/2 



1 + 7^(fc + qCOS0Kq) 



Kf - . 

hp 

w k + qcos6Kq (48) 
which is valid to leading order in A/Kp, and where 

cos6lA'g ^ k-q (49) 

^ r cosiOK - 0q) 

y cos 9k cos 9q + sin 6k sin 9q cos{ipK ~ Vq) 

(50) 

in d = 2 and c? = 3, respectively. Equation P7)) now 
becomes 



bf 



z + 2-d 

s 9 



d'^q'dk'd'^-^Vljfde'duu'i^'ip'-^' 



x9( A/s — \k + qcos9Kq\ 
A 

d'^q'dk'd''-^nr>de'duj'xj/tl:'(f>' 



z + 2-d 

s 2^ g 



xe(^A- \k' + s^'-^'^^'q' C0s9Kq\j 



(51) 
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where k' = sk, q' ~ s^^^q, e' = se, and uj' = slo. Clearly, 
for generic values of z the 0-function does not return to 
its original form after the renormalization group trans- 
formation. We should be pleased, however, that in the 
special case z = 1, the 0-function is form-invariant. 

e(A/s-|/C2|) « eiA/s~\k + qcoseKq\) 

= e{A- s\k + qcoseKq\) 

= e{A - \k' + q' COS 9Kq\) (52) 

The boson-fermion coupling can now be written 

f-A 



s: 



bf 



3-d 

S 2 g 



cfiq'dk'd'^-^VLgde'duj'iP'iP'cj)' 



xe A - \k' + q' coseKq\ 



and we can identify 
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,(3-d)/2 
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which is equivalent to 



[g] = (3-d)/2 



(53) 



(54) 



(55) 



This is one of the central results of this paper. The 
coupling is marginal in d = 3 and relevant in d = 2. 
Of course, this result depends on the choice of field di- 
mensions; equations (|46p with z = 1. In application 
to an antiferromagnetic Kondo lattice, we have previ- 
ously developed a model where the boson dimension is 
—d [rather than equation PS|) ] and used the scheme ex- 
plained here to show that the boson-fermion coupling is 
exactly marginal in that case^^. 

Equation (|55p is only valid when z = 1 because only 
then is the 0-function form-invariant. What can be done 
when z 7^ 1? This question is particularly pertinent 
to the controversy surrounding the renormalization of a 
gauge-field coupled to a fermion with a Fermi surface 
^ji0jiAji2_ It is also germane to the ferromagnetic phase 
of heavy fermion systemsSi. Section IVII will present a 
different scheme that is applicable to problems with arbi- 
trary values of z and also reproduces equation (|55p when 
z = 1. Here, we merely explain why the present scheme 
fails when z ^1. 

We have actually already seen the problem in equation 
(|5ip . where it is obvious that the 0-function is not form 
invariant. This is similar to the dilemma we encountered 
in the pure fermion problem, as seen in equation (|32p . To 
make progress, we try the same strategy used in the pure- 
fermion problem where we restricted our consideration to 
the phase space where the 0-function does scale perfectly. 
We did so by demanding /C4 = (or | A| = 1) in ^4 , which 
can be implemented by simply inserting 5{K,4). Here, the 
analogue of that additional constraint is IC2 = 0. This 
new condition can also be written 



Thus, besides staying within their respective cutoffs, the 
choices available to K and q , when IC2 = 0, are restricted 
in such a way that their sum vector must sit precisely 
on the Fermi surface. Once K is chosen, q is obligated 
to connect K + q to the Fermi surface thus limiting its 
permissible magnitudes and angles quite severely. This 
is depicted in figure [21 




FIG. 2: K must stay within the annulus while q must stay 
inside the little circle of radius A. Under the restriction IC2 = 
0, the sum K + q must sit precisely on the Fermi surface. 
The only phase space that satisfies /C2 = is the thick gray 
line which represents a small patch on the Fermi surface of 
size 0{A'^~^). Clearly, the limit /C2 = is not the same as 
A/Kf ~ since the latter would shrink the gray patch to a 
point. 

Under the restriction IC2 = 0, the boson-fermion cou- 
pling can be written: 

S^/ = s^^g / d'^q'dk'd'^-^nj^de'duj'ip'ip'cf)' 

xO(A/s- |/C2|)(5(/C2) (57) 

2-d 

= s^^g / d'^q'dk'd'^-^Vlj^deduj'ij'ij'cl)' 



xj(fc' + s(^-i)/^9' cos6Ik,) 



\K + q\ = 



(56) 



(58) 



This should be compared with equation ([39)1 . In the pure 
fermion case, we showed there that whether or not we in- 
sert 5[K,4) makes no difference to the value of [u f] because 
(5(/C4) is ^ (5(6*3 — 6*1) or ^ S {63 — 62), which has zero scal- 
ing dimension. Furthermore, this additional constraint is 
of a non-singular nature. 

In contrast, for the boson-fermion coupling in equation 
(|58p . the insertion of (5(/C2) involves a dimensionful quan- 
tity. If we were to integrate against (5(/C2) and eliminate 
6q as suggested by figure [U we would induce an addi- 
tional 1/momentum factor in violation of the RG edict 
that the coupling be a non-singular function of momen- 
tum. More intuitively, figure [5] shows that imposing the 
constraint IC2 = singles out an unrealistic sort of cou- 
pling that glues the out-going fermion iI'k^^ to Fermi 
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surface regardless of the value of K or q. This no longer 
represents a generic forward scattering process, and is of 
no interest to us. How to correctly capture a generic for- 
warding scattering process will be discussed in section lVTl 
At this point, a few issues are worth emphasizing. 

• Since we integrated against the delta functions in 
equation (j4T|) . energy and momentum are explicitly 
conserved. 

• The quantity IC2 is not a free variable and it does 
not necessarily scale in the same way as bosonic 
or fermionic momenta. This is consistent with 
the non-scaling of /C4 in the pure-fermion problem 
when |A| 7^ 1. Only when z — 1 docs JC2, and thus 
the constraint 0(A — |/C2|)i scale in a simple way. 

• In this scheme, all components of q scale the same 
way. In particular, [df^q] = d/ z = d. At the same 
time, only fermionic momenta in the direction nor- 
mal to the Fermi surface scale. 

• Here, k is not a vector. It does not have parallel 
or perpendicular components as discussed in cer- 
tain patching schemes. For more on the patching 
scheme, see section IVll 



• Although the RG scheme developed in this section 
does not work for general values of z, it is perfectly 
well suited to the special case z = 1. 

• Figure [2] gives us an important hint about what 
may be happening for z > 1. Since k' ~ sk and 
q' = s^^^q, we know that after several iterations 
of the RG, the deviation of K from the Fermi sur- 
face will be much smaller than the magnitude of q, 
i.e. k <^ q. As a result q will tend to point in a 
direction perpendicular to K, which means it will 
be very nearly tangent to the Fermi surface. In 
this way, it may seem as if bosonic momenta scale 
anisotropically in a local coordinate system defined 
with respect to the direction determined by a fixed 
K. This important observation will be developed 
more fully in section IVTl when we devise a scheme 
suitable to z 7^ 1. 



V. CHOICE OF MOMENTUM INTEGRATION 

Before moving on to the general case z 7^ 1 , in this sec- 
tion we resolve a seeming ambiguity for the scheme we 
developed in the previous section. As we found in equa- 
tions (|^ and there are two ways to express in 
momentum space. We have already shown in detail that 
making the choice in equation can yield a consistent 
RG prescription. Now we will show why the alternative 
decomposition 



is not an appropriate starting point to determine the scal- 
ing dimension of the boson-fermion coupling. The prob- 
lem is that the argument of the boson field, Q = A'l — A'2, 
does not transform homogeneously, so we do not know 
what dimension to assign to the boson itself. To see 
this, write each fermion momentum vector in terms of 
a direction and a deviation from the Fermi surface: 
= [Kp + ki)ki. This gives 



ISI = 



[Kl + Kl - 2K1K2 cos 612] ^'"^ 

(Kf + kif + {Kf + k2f - 2{Kf + kl) 
1 1/2 

x{Kf + k2) cos 6*12 



KfV2 



fl — cos ( 



kl + k2 



1/2 



« KF^/2il-cos9i2) (1 + ^^^) (59) 
which is true to leading order in 1/Na, and where 

cos 6*12 = Kl ■ K2 

cos(0i - 62) 

cos 9i cos 62 + sin 9i sin 9i cos{ipi — ^2) 

(60) 

for d = 2 and d = 3. In equation ([23]) we committed to a 
specific prescription in making scale invariant where 
angular components of the momentum do not scale. We 
therefore cannot allow angles to scale in cither. Using 
the specific prescription in equations (|46p determined by 
the quadratic parts of the action we find 

S\'^ = f s^^dk[s^^dk2S^^de'iS^^de2 



s^^^tp' {Kf + k'2,i€'2)s^^^tlj'{KF + k[,ie'i) 



X(/.fA:FV2(l-cos6li2) 



1 



2sKf 

s ^z4-zs-M)5(2,l)e(AAsi/^-|Q|) 



(61) 



where 



e(A/si/^-|Q|) = e A-.si/^ifj.V2(l-cos0i2) 



k-^ ~\~ 
2sKf 



(62) 



d'^Kid'^K2 



^K^^K^'I^K,-K,9{K2,Ki)Q{K^\Q\) 



Notice that in equation (|6ip the fermion fields arc primed 
whereas the boson field is not. 

There are two problems. First, the 8-function does 
not return to its original form, making it impossible to 
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compare the flow of tlie coupling function before and after 
the RG transformation. This is the same problem we 
encountered in the pure fermion case of section [TTTl and 
the other boson- fermion prescription from section [TVl 

Second, we have a new dilemma, which is that we do 
not know how the 4> field transforms under the change of 
argument in equation (j6ip . All we know from equation 
(fTTjl is that 



(63) 



which states that the boson scales in a (generalized) ho- 
mogeneous fashion. If we transform the boson arguments 
in a non- homogeneous way. as in (jHT]) . we are not guar- 
anteed that such a coordinate transformation will induce 
a simple multiplicative prefactor. Note that the mathe- 
matical requirement that the boson field transform ho- 
mogeneously means that the relative angle between the 
incoming and outgoing fermions must be allowed to scale. 
Said another way, when the magnitude of q scales, the an- 
gle of K2 = {Ki + q) must change. However, when we 
choose to work in representation (j43p . all momenta are 
fermionic which forces the wrong type of rescaling on the 
boson field. 

Thus, we cannot adopt representation (|43p because, 
first, the 0-function is not form-invariant, and second, it 
forces a non-homogeneous coordinate transformation on 
the boson field. 

How might we try to remedy these two problems? We 
could attempt the same strategy that worked in the pure 
fermion case where we restricted our consideration to 
/C4 = 0; see equations ((33)) - ([35|) . However, when both 
bosons and fermions are present this tactic is bound to 
fail. We already found this in section HVl where we con- 
sidered the limit /C2 = using representation (|32|). Here, 
the analogous restriction is Q = 0. Under these circum- 
stances, equivalent to Ki = A'2, the 0-function is triv- 
ially invariant. However, the boson loses its field charac- 
ter with 0(0) not scaling at all. 

Let us attempt a different remedy by relaxing the re- 
striction slightly and consider Ki = K2. This is equiva- 
lent to Q = \ki — k2\- Here, the 0-function will not be 
form invariant because it transforms to 



e(A 



\ki^k2\) = e(A 



Undaunted, we make the further restriction to z = 1, 
in which case the 0-function is form- invariant, and now 
less trivially so. Unfortunately, we still have the problem 
that the boson field scales unnaturally: (j){\Q\,6q,Lpq) ~ 
(j){\ki — fc2 1, const, const). As show in figure [31 the boson 
momentum vector, which is defined as the vector joining 
the tips of Ki and K2, lies directly parallel to Ki = K2- 
When we force Ki = K2, the boson momentum loses 
its angular freedom and thus no longer scales homoge- 
neously. 

One might wonder why we are being so strict about 
the form of the field transformation when it seems like 




FIG. 3: Under the restriction Ki = K2, the boson momentum 
q can still have a non-zero modulus, but it loses its angular 
freedom and is forced to point exactly normal to the Fermi 
surface. This constraint does not allow the boson field to 
transform in a homogeneous fashion. 



the other scheme in equation (|^^ 



d^Kd^q ^g^.ijj^cj)q^g{K, q)e{A - I/C2I) (64) 



also violates this principle. In fact, the fermion is not re- 
quired to be a homogeneous function of momentum any- 
way. All that we need from equation ([23)) is 



^{Kp + s-^k') = s^/^^'{Kf + k') 



(65) 



The incoming fermion in equation (|64p is clearly of this 
form, whereas the outgoing fermion can be written: 

i;{\K + q\) « tlj{KF + s-^k' + s-\' cosOKg) (66) 
In this form, we know this expression is equivalent to: 



ll^iKp + S" fc' + S~^q' COS0Kq) 

s^/'^iP'{Kf + k' + q' cosOKq) 



(67) 



Thus, both fermions in equation ()64p transform as equa- 
tion ([^5)) . Finally, (f){q,iuj) obviously scales according to 
equation ([63)) [and ([TT)) ]. Therefore, in equation ([64)) we 
know how all fields transform under equations (|46p . so 
representation ()42p suffers from none of the shortcom- 
ings we identified for representation ([^5)) . 

With this new understanding, we should also check 
that in the pure fermion problem the field ipilC^) trans- 
forms in a consistent manner. To see this, we need to 
keep a few more higher order terms than what we showed 
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earlier. 

|i?4|2 = + Kl + Kl + 2Ki ■ K2 - 2i?i • K3 

-2K2 ■ K3 
w 2KF{ki +k2 + h + SKpfi) 

-2KfKi ■ k^iki +k3+ Kf) 
-2A>A'2 • A'3(fc2 + /ca + A» 
= 2KF[KF{ki ■ A'2 - ATi • A'a - A'2 • A'3 + 3/2) 

+fci[i + A'i-(A'2-A:3)] 

+fc2[l + A'2- (A'l- ATa)] 
+fc3[l-A:3- (-ft:i+A'2)]} 
= 2Kf\^Kf\M^ /2 

+ki[i + kl ■ {k2 - ks)] 
+k2[i + k2 ■ {kl ~ ks)] 

+k3[l ~ ka ■ {kl + k2)]} (68) 

where we used ([31]) . In the special case where ATi = ATa, 
corresponding to forward scattering, we have 

\ki\ « Kf + k2 + {kl - A:3)A'i • A:2 (69) 
This shows that 

^(|i?4|) = ^{KF + S-'^k'2+S-\k'i-k'3)ki-k2) 

(70) 

which is precisely the scaling form appropriate for a 
fermion in equation (pS)) . In the same way, it is easy to 
show that the fermion scales appropriately for the cases 
k2 = A^3 and ATi = — Ar2. In these cases we have 

|i?4| ^ Kf + kl + {k2 - k3)ki ■ ks (71) 
I .^4 1 « KF + k3 + {k2- ki)ki ■ k3 (72) 

respectively. Thus, all the results of Shankar remain 
valid. 

Finally, it may at first seem puzzling that the scal- 
ing of the constraint in equation (|62p is so problematic 
since we were able to find a simple solution in the pure 
fermion problem involving S[. There, the constraint in- 
volved IC4 = \kz — k2 — ki\ — Kf, which measures a 
deviation from the Fermi surface. However, in the repre- 
sentation of the boson- fermion coupling in equation P5)) . 
the constraint involves Q = A'2 — ATi which is not a de- 
viation from the Fermi surface and as written, can take 
any value between and 2Kf\ see equation ((55)) . 

To summarize, '0(^4) scales like a fermion, ^{K,2) 
scales like a fermion, but (^{Q) does not scale like a boson. 
We therefore cannot use equation ([43]) to represent the 
boson-fermion coupling because we do not have knowl- 
edge of the boson field scaling under such a coordinate 
transformation. 



VI. THE PATCHING SCHEME 

When 2; 7^ 1 , the scheme we developed in section IIVI 
no longer works. The problem is that under mode elimi- 
nation and rescaling, the constraint function changes its 
form. 

e(A/s-|/C2|) « Q{Kls-\k + qcosBKci\) 
= Q{K- s\k + q cose Kq\) 
^ e(A- |fc'-f s^^'i^/^g'cos^if,!) 

(73) 

where k' = sk and q' = s^l^q. When z 7^ 1, we cannot 
reliably determine the flow of the coupling because the 
structure of the interaction itself has changed under this 
RG transformation. This is the same dilemma encoun- 
tered in the pure fermion problem in equation ([5^ . Also, 
notice that writing 53'' in terms of an integral over Ki 
and A'2, rather than AT and q, will not cure the problem. 
In the previous section we explained why this is the case, 
even for z = \. 

For these reasons, when 2; 7^ 1 we adopt a different 
method where we scale toward a specific point on the 
Fermi surface. Although the details differ, this is similar 
in spirit to some previous work on the renormalization of 
the gauge-spinon problem^ '^^i^^'^^ . 

In d = 2, consider the annulus in momentum space 
defined by — A < A: < A. Now subdivide the annulus 
into A'a regions of angular size A0: N\A9 = 2tt 
A9 = 2ttA/Kf- Each patch will be approximately of size 
^ A^. The same idea is easily generalized to c? > 2. This 
should be familiar from multidimensional bosonization^ 
and functional R G^^i^^i^^'^^ ; we refer the reader to those 
papers for further details. 

The momentum integral for the quadratic part of the 
fermionic action, S2 , can now be decomposed into a sum 
over A'a identical patches: 

si = / d'^Kdei;{ie-^g)i; 

Na „A 

= J2 d'^kpdepil}p{iep-£,^^)i}p (74) 

Here, kp = (fcp__L,fcp ||) is a local coordinate within each 
patch which has components parallel and perpendicu- 
lar to some reference frame. We define this special lo- 
cal reference direction to be the normal vector to the 
Fermi surface at the patch origin. Thus, is tangent 
to the Fermi surface at the patch origin. Within each 
patch, functions of momentum can be expanded around 
the patch origin. Consider, for example, the patch cen- 
tered at AT = (0, Kf) which we will label as patch p = I, 
and where we have specialized to d = 2 for concreteness. 
Near this point, the dispersion of a perfectly parabolic 
band can be expressed in terms of local patch coordi- 
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nates as follows. 



■{O.Kf) 



{K 



.'2 



+ ^'y)l(0,ifF) 



1 

2m 



2m 2 



2m 
2Kf 



0) 

iff) 
-0)2 



2m 



if,, 



(75) 
Kp^ and 



where for this particular patch, fc^, 
fcij^ = Kx- We have also defined vp = Kp/m and a = 
vp/{2Kp) ~ l/(2m). As a sum over all the patches that 
enclose the Fermi surface, the quadratic part of the action 
can now be written 



S. 



f 



Na 

E 

p=i 



i'^kpdepipp (^iep - vpkp^\\ - afc^ ipp 



Note that the concepts of parallel and perpendicular only 
make sense with respect to a perfectly flat surface, or the 
normal to a specific point on a curved surface. Wc take 
this specific point to be the center of the patch. Mo- 
mentum components in the same direction as the vec- 
tor normal to the Fermi surface at the patch origin are 
considered "parallel," whereas momenta tangent to the 
Fermi surface are labeled "perpendicular." We caution 
that different conventions exist in the literature for what 
is deemed parallel or perpendincular. We adopt the con- 
vention of referenceii. 

Within each patch, the momentum integral is limited 
to a box of dimension A in every direction. For example, 
in d = 2 this means 



(fk = 



dk 



-A 



dk\ 



-A 



(76) 



We have dropped the patch indices since wc assume all 
patches are identical in the sense that variables scale in 
the same manner in every patch. In the absence of van 
Hove singularities and nesting instabilities, this is a rea- 
sonable assumption—. 

Within this patching formalism, and when we consider 
only (f « so that the entire boson phase space can be 
restricted to a single patch, the quadratic part of the 
bosonic action can be written in straightforward fashion. 
To be concrete, consider z = 3: 



d'^qduj 




Within each patch, bosonic momenta q and fermionic 
momenta k are all measured with respect to the same 
single point, the patch origin. Consequently, bosonic and 
fermionic momenta scale the same way, that is 



(77) 
(78) 



See figure 21 Whether wc label momenta by fc or g is 



k + q 




FIG. 4: Local coordinate system in patch p = 1 whose center 
is located aX K = (0,AV). Bosonic and fermionic momenta 
now scale identically 

thus immaterial since they scale identically; this is in 
stark contrast to the scheme developed in section Hvl for 
the case z = 1. 

Of course, the possibility exists that [fcy] ^ [k±]- In 
fact, we will now argue why they cannot be the same 
when z 7^ 1. 

The fixed point is defined by constructing the scaling 
scheme so that the quadratic part of the action, S2+ S^, 
is scale invariant. Scale invariance of requires 

[e] = K] + [fc||] = [a]+2[fcj 



(79) 



while scale invariance of S2 necessitates 



[7] 



z/2' 



(80) 



Next, we observe that since we want to scale bosons and 
fermions simultaneously, it is sensible to give them equal 
scaling dimensions which wc denote by 



[E] 



(81) 



The conditions from the quadratic parts of the action 
now become 



[E] 



[E] = [vf] 



-2 I 2 
<1± + 111 



[fc||] = [a] 

z/2 



(82) 
(83) 
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At this point we demand that the dispersion relations of 
low-energy excitations be preserved under scaling. For 
fermions near the Fermi surface, energy must be a linear 
function of momentum. We thus set [vp] = to establish 



[E] = [fc||] = [a]+2[fcj 
This is furthermore justified by the fact that 



(84) 



Vpka 



< 1 provided fc^ji < A < Kp. The "curvature" 

term is a small correction. Thus, parallel momenta scale 
like energy. 

We now use ([84]) and ((78|) in equation ([83]) to determine 
the dimension of the perpendicular momentum: 



z/2 



In the infrared limit, this becomes 



(85) 



(86) 



because E^^^ > E^ when z > 1. To preserve the bosonic 
dispersion [i.e. w ~ q^), we set [7] = 0, obtaining 



[E]/z 



Now we plug this result into equation 
dimension of a: 



(87) 
to find the 



[E] 



[a\ 
(1 



- 2 [E] /z 
21 z) {E\ 



Finally, we are free to choose the value of \E\, which 
we set equal to unity for convenience; any other value 
will only induce the same multiplicative prefactor on all 
dimensions, but relative dimensions will be unaffected. 
To summarize. 



M 
k 

[a\ 
[vf] 
[7] 



[e] = [c] = 1 
[911] =[E] = l 
[gl] - [E] /z = 1/z 

1 - 2/z 







(89) 



Note that for 2 = 3 the dimension of fcj^ appears to 

suggest that the fermionic bandstructure changes under 
scaling. This is an illusion since the parameter a, which 
is a measure of the curvature, is allowed to scale in order 
to precisely compensate the scaling of k± , thus ensuring 
that the band remains invariant. 

Plugging these values into the quadratic action yields 
the dimensions of the fields: 

m = -'-^^ (90) 



2z 



2z 



1 



2z 



(91) 



We now have enough information to determine the di- 
mension of the boson-fcrmion coupling. 



s: 



bf 



J (fkid'^k2d'^qde2deidw 

patches 

xg{k2, ki,q, 62, ei,uj) i/i(2)V'(l)^(g, oj) 
y-S^^-^Khi- - ki.± - q±)S{k2,\\ - fci.ii - <7||) 
X(5(e2 - ei - t^) 
xe(A-|fc2,±|)e(A-|fc2,|||) 

xe(A-|fci,i|)e(A-|fci,|||) 

xe(A-|gl|)e(A-|g|||) (92) 

Note that this is slightly less general than could be the 
case. We have restricted our consideration to nearly for- 
ward scattering processes which means that (f w or, 
equivalently, ki and k2 belong to the same patch. Inter- 
patch processes, such as the BCS instability, are not cap- 
tured. 

Since we are scaling toward a single point, momentum 
and energy conserving delta functions and cutoff con- 
straints factorize nicely. Integrating against the delta 
functions yields 



J d'^kd'^qdeduj 

patches 

xg{k, q, e, uj) Tp{k + q,e + uj)%p{k, e)(j){q, uj) 
xe(A- |fci + gl|)e(A- +g|||) (93) 



where we have placed some of the constraints in the limits 
of integration. Unlike what happened in section llVl there 
is no difference in eliminating boson or fermionic vari- 
ables due to equations (|89p . Here it is arbitrary whether 
we call momentum k ov q since in the patching scheme 
they scale the same way. Additionally, the factorization 
of parallel and perpendicular components of momenta 
means the arguments of the fields scale in a straightfor- 
ward fashion. Indeed, after mode elimination and rescal- 
ing we find: 



bf 



9 E 

patches 



xs^ 



{■iz+d-l)/{2z)^l^{2z+d+l)/{2z)^l 



xe(A/,sl/- - ,s-l/-|fc^ + 9'^|)e(A/s - s-'\k[ + q[\) 



„(3-<i)/(2 



d'^k'd'^q'de'duj'^'^'cj)' 



patches ' 

xe(A-|fc^ + gll)e(A-|fcf|+9[||) 



(94) 



where we have Taylor expanded g and kept the most 
relevant (constant) piece. In this patching scheme, the 
constraints and fields transform in a simple way, so we 
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can simply read off the dimension of the eouphng. 



[<?] = 



3-d 
2z 



(95) 



The relevance or irrelevance of this coupling is in some 
sense arbitrary outside the context of a specific physical 
problem. The value of [g] depends crucially on the di- 
mensions [</)] and ["0] , and these will be determined by the 
problem under consideration. For example, in the con- 
text of magnetic phases of the Kondo lattice, se o^^i^^'^^ . 
Several important comments are now in order: 

• The result in equation ([M)) is identical to when 
z = 1. Therefore, the patching scheme developed in 
this section yields an answer equivalent to the ex- 
tension of Shankar's scheme presented in section llVl 
using global coordinates. While the latter approach 
is perhaps more intuitive, it is not justifiable when 
z 1. On the other hand, the patching scheme 
requires same careful interpretation, as discussed 
below, but is consistent for any value of z. 

• It is necessary to give the curvature parameter, a, a 
non-zero scaling dimension in order to compensate 
for the way that k± scales. Rest assured, however, 
that [afc^] = [wf^||] = [e] so that the fermion band 
is kept invariant. In this way, we do not need to 
scale the number of patches. 

• It may seem as if the bosons have become 
anisotropic, but this is an illusion due to the nature 
of the local coordinates we have chosen. Because 
of the sum over patches, we have included an equal 
weighting of q components in all directions, even 
though locally we only keep q± within each patch. 
Of course, it does mean that in the low-energy limit 
bosons become locally tangent to the Fermi surface 
for fixed value of fermionic momentum K . This is 
not surprising and was noticed long ago^^. We 
even saw hints of this in section HVl In that scaling 
scheme k' = sk and q' ~ a^l^q. When z > 1, the 
length of \q\ scales more slowly than the deviation 
from the Fermi surface, k. As a result, in the low 
energy limit, the boson momentum will tend to lie 
tangent to the Fermi surface. 

• In the patching formalism, the dimension of the 
boson field in equation (j9ip derives from 



and similarly for the fermion field. Once again this 
takes the form of a generalized homogeneous func- 
tion, but is different from the type of scaling in 
equation pT|) or (|23p . For more on generalized ho- 
mogeneous functions, see^^. 

• The form of the interaction we consider is limited 
to nearly forward scattering (g k, 0) intra-patch 
processes. Inter-patch processes are not captured. 



and this makes comparisons with the pure-fermion 
RG somewhat delicate. Consider a four-fermion in- 
teraction with incoming momenta K\ and K^-, and 
outgoing momenta Kj, and The difference be- 
tween incoming and outgoing momenta at the left 
vertex can be small, say — K\ = qidt vertex ~ 0. 
This can match up with small momentum transfer 
on the right: K^- K2 = bright vertex ~ 0. However, 
this says nothing about the relationship between 
Ki and K2- Indeed, Ki and K2 can each inde- 
pendently take any value around the Fermi sur- 
face, i.e. \K2 — Ki\ can take any value between 
and 2Kf. Thus, "forward scattering" processes 
in a boson-fermion formalism are not necessarily 
equivalent to "forward scattering" processes in a 
four-fermion formalism. The latter (four-fermion 
coupling) involves two patches, whereas the former 
(boson-fermion coupling) involves only one patch. 
In other words, the dimension of u/ is not simply 
given by [g^]. 

• If we were to include self energy corrections into 
Sl and establish this as the new fixed point, the 
values of the dimension assignments would change, 
but the philosophy would be the same. For ex- 
ample, in the gauge-spinonii and ferromagnetic 
Kondo lattice systems^!, gapless overdamped z = ?> 
bosons lead to a characteristic electron self energy 
I](e) - £2/3 in d = 2 and S](e) - -eloge in 
d = 3. We can define the new fixed-point action 
with S2 = j %lj{e'^/^ — wf^ii — ak'^)ip. Using the 
same philosophy defined in this section, we would 
assign 

[E] = [e] = [u] = z/d 



[a\ 
[vf] 
[7] 



= [E] /z = 1/d 

1 - 2/d 







(96) 



This also leads to a change in the dimensions of 
the fields and the couplings, but the methodology 
is no different than what has already been discussed 
above. See^i for further discussion. 

• There is some debate in the literature about how to 
properly scale the gauge-spinon model which corre- 
sponds to z = 3 91IO1I1, 12^29 ^ rpj_^^ consistent scaling 
scheme within the patching formalism we advocate 
here coincides with that ofii. 



VII. CONCLUSION 

This paper has developed an easy-to-use RG proce- 
dure for theories containing both bosons and fermions 
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with a Fermi surface. We reviewed the global coordi- 
nate approach to the fermionic RG as formulated by 
Shankar, showed how to generalize this formalism to in- 
clude bosons with dynamical exponent z = 1, and ex- 
plained why such an approach will not work when z ^ 1. 
We pointed out that a consistent scheme must ensure 
that the kinematic constraints, which result from the 
conservation of momentum and the effective field the- 
ory cutoffs, remain invariant to the RG transformation. 
In addition, field rescaling can only be properly identi- 
fied in interaction terms when the coordinates of the field 
transform in a known way, as specified by the quadratic 
part of the action. 

We also showed that, for z = 1, the same results arise 
within a patching scheme. Here the momentum space 
near the Fermi surface is partitioned into patches. For 
z 7^ 1, the patching scheme represents the only consistent 
RG approach to mixed fcrmion-boson systems. 

Coupled boson and fermion problems arise in a variety 



of contexts. We have already mentioned the problems 
of itinerant magnets which have directly motivated our 
work here, as well as the subject of gauge fields coupled 
to fermions. In addition, fermion-boson mixtures of cold 
atomic gases^ may provide another interesting setting 
for this work. We hope the RG program described here 
will be useful for related problems in other settings as 
well. 
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